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HIGHLIGHTS 


►  Analysis  of  the  combined  effect  of  multiple  degradation  processes  in  SOFC  stacks. 

►  Acceleration  of  the  degradation  and  sequential  activation  of  processes  is  predicted. 

►  Optimisation  of  the  operating  conditions  can  extend  lifetime  by  a  factor  up  to  10. 

►  Operating  conditions  for  the  highest  performance  and  the  best  durability  differ. 

►  Overpotential  rather  than  current  density  influences  the  degradation  behaviour. 
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The  degradation  of  solid  oxide  fuel  cells  (SOFC)  depends  on  stack  and  system  design  and  operation.  A 
methodology  to  evaluate  synergistically  these  aspects  to  achieve  the  lowest  production  cost  of  electricity 
has  not  yet  been  developed. 

A  repeating  unit  model,  with  as  degradation  processes  the  decrease  in  ionic  conductivity  of  the 
electrolyte,  metallic  interconnect  corrosion,  anode  nickel  particles  coarsening  and  cathode  chromium 
contamination,  is  used  to  investigate  the  impact  of  the  operating  conditions  on  the  lifetime  of  an  SOFC 
system.  It  predicts  acceleration  of  the  degradation  due  to  the  sequential  activation  of  multiple  processes. 
The  requirements  for  the  highest  system  efficiency  at  start  and  at  long-term  differ.  Among  the  selected 
degradation  processes,  those  on  the  cathode  side  here  dominate.  Simulations  suggest  that  operation  at 
lower  system  specific  power  and  higher  stack  temperature  can  extend  the  lifetime  by  a  factor  up  to  10, 
because  the  beneficial  decrease  in  cathode  overpotential  prevails  over  the  higher  release  of  volatile 
chromium  species,  faster  metallic  interconnect  corrosion  and  higher  thermodynamic  risks  of  zirconate 
formation,  for  maximum  SRU  temperature  below  1150  K.  The  counter-flow  configuration,  combined  with 
the  beneficial  effect  of  internal  reforming  on  lowering  the  parasitic  air  blower  consumption,  similarly 
yields  longer  lifetime  than  co-flow. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Distributed  power  generation  is  a  foreseen  application  for  solid 
oxide  fuel  cell  (SOFC)  devices.  Combined  heat  and  power  systems, 
even  small  ones,  based  on  this  technology  could  achieve  high 
efficiencies,  fuel  flexibility  and  load  following  capabilities  along 
with  low  maintenance  costs  and  noise  nuisance.  Despite  many 
advances,  SOFC  prototypes  do  not  yet  meet  the  reliability  and 
durability  standards  of  the  industry.  Unresolved  structural  issues 
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affect  the  reliability  by  hindering  a  completely  safe  assembly  and 
operation  of  the  cells,  whereas  the  degradation  of  the  electro¬ 
chemical  performance  limits  the  durability. 

SOFC  system  design  and  operation  must  aim  at  the  lowest 
achievable  cost  of  electricity  under  practical  conditions  that 
support  utility  grid  operation,  as  well  as  end-user  and  owner 
demands  [1].  To  appropriately  meet  these  requirements,  they  have 
to  be  integrated  with  lifetime  analysis,  since  they  are  known  to 
have  a  strong  impact  on  the  SOFC  degradation  [2,3].  Because  of  the 
various  underlying  physical  phenomena,  this  involves  to  retrieve 
knowledge  from  different  research  fields  and  length  scales.  To  our 
knowledge,  an  integrated  methodology  to  evaluate  stack  and 
system  design  and  operation  in  the  view  of  electrochemical  and 
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mechanical  degradation  has  not  yet  been  developed.  Such  a  tool 
could  help  to  further  approach  the  potential  capabilities  of  the  SOFC 
technology  [4,5]. 

Button  cell  endurance  tests  in  inert  alumina  test  arrangements 
exhibit  degradation,  induced  by  localised  microstructural  alter¬ 
ations  of  the  materials  and  interfaces  of  the  membrane  electrode 
assembly  (MEA),  due  to  particle  coarsening  [6-11],  detrimental 
reactions  between  the  materials  [12-14]  or  impurities  in  the 
starting  powders  [15].  In  an  SOFC  stack,  the  cost-effective  inter¬ 
connects,  sealants  and  gas  diffusion  layers  (GDL)  amplify  material 
incompatibilities,  hence  the  degradation,  because  of  their  lower 
stability  than  the  MEA  materials  in  the  SOFC  environment  [16].  The 
transport  of  contaminants  by  the  gases  increases  the  distance  over 
which  detrimental  interactions  operate  [17].  The  uneven  distribu¬ 
tions  of  temperature,  current  density,  gas  composition  and  over¬ 
potential  in  a  stack,  induce  different  local  degradation  rates.  The 
interactions  between  the  detrimental  processes  govern  their 
evolutions,  which  accounts  for  a  progressive  acceleration  of  the 
degradation  rate  [18].  An  efficient  mitigation  strategy  must  deter¬ 
mine  a  suitable  nominal  operating  point  that  balances  power 
density  and  efficiency  requirements  with  lifetime.  Spatial  temper¬ 
ature  control  can  alleviate  the  progressive  accumulation  of  damage 
under  comparatively  more  harmful  operation  during  load 
following,  which  can  affect  the  durability  of  the  SOFC  stack  [19]. 

Owing  to  the  multiphysics  nature  of  degradation,  models  are  of 
interest  to  first  improve  the  relevance  and  efficiency  of  the  develop¬ 
ment  cycles,  and  then  to  adjust  the  product  for  each  application.  The 
mitigation  of  degradation  through  improvements  of  the  stability  of 
the  MEA  and  other  single  repeating  unit  (SRU)  components  is  in 
progress  and  controversy  remains.  As  discussed  in  Ref.  [20],  because 
the  focus  of  experimental  studies  is  principally  on  either  (i)  gaining 
knowledge  on  the  underlying  processes  or  (ii)  scanning  for  the  best 
technological  solution,  the  range  of  investigated  conditions  does 
usually  not  cover  that  found  in  SOFC  stacks  in  operation.  Flerein  lies 
their  partial  incompatibility,  in  particular  that  of  (ii),  with  modelling 
at  the  SRU/stack  scale.  The  range  of  operation  is  determined  by 
indicators  for  anode  reoxidation  [21  ]  or  cathode  reduction  in  the  case 
of  imperfect  sealants  [22,23].  The  SRU  geometry  and  gas  manifold  is 
adjusted  to  reduce  high  local  current  densities  to  avoid  ensuing  high 
localised  rates  of  heat  generation  and  current  constriction  effects. 
More  refined  analyses  [24]  seek  to  identify  the  patterns  of  degrada¬ 
tion  in  electrochemical  impedance  spectra  (EIS)  for  in-operation 
diagnostics.  The  next  step  is  to  implement  the  multi-factorial  and 
history-dependent  nature  of  SOFC  degradation  in  SRU  models,  for 
reliable  lifetime  predictions  and  integrated  mitigation  approaches. 
The  state-of-the-art  understanding  of  the  SOFC  degradation  is  not 
precise  enough  to  support  the  completion  of  this  difficult  aim.  It 
however  provides  a  starting  point  that  can  serve,  in  synergy  with 
experimental  tasks,  to  gain  knowledge  on  the  underlying  physics  and 
must  be  accordingly  and  continuously  updated. 

In  this  study,  a  simulation  tool  has  been  developed  to  support 
stack  design  and  operation  in  the  view  of  both  lifetime  and  elec¬ 
trical  efficiency.  Our  calibrated  electrochemical  model,  which 
includes  several  simultaneous  degradation  phenomena  [20,25]  is 
implemented  in  SRU  models  of  a  planar,  intermediate-temperature, 
anode-supported  stack  design  to  investigate  their  combined 
detrimental  effect  on  the  durability  under  the  practical  operating 
conditions  an  SOFC  stack  must  withstand  in  the  field.  The  decrease 
in  ionic  conductivity  of  8YSZ,  metallic  interconnect  (MIC)  corrosion, 
nickel  particle  growth  and  chromium  contamination  of  the 
LSM-YSZ  cathode,  completed  by  transport  of  volatile  chromium 
species  [22],  are  considered.  This  Part  I  investigates  the  effects  of 
operating  conditions  on  the  start-up  and  long-term  performance  of 
SOFC  stacks,  to  identify  the  potential  of  this  practical  approach  to 
extend  the  lifetime  and  assist  the  exploration  of  the  operating 


envelope  of  SOFC  devices.  Part  II  [26]  focuses  on  the  distribution  of 
the  degradation  over  the  active  area.  System  considerations  are 
restricted  to  the  parasitic  loss  of  the  air  blower  and  efficiency  of  the 
inverter,  to  avoid  the  computational  burden  of  a  complete  inte¬ 
grated  system  model,  yet  capture  the  key  system  efficiency  char¬ 
acteristics.  Operating  conditions  are  optimised  for  co-flow  and 
counter-flow  configurations,  along  the  trade-off  between  system 
power  and  efficiency,  first  for  the  best  initial,  then  for  the  best  long¬ 
term  electrical  efficiency  during  operation  at  a  fixed  system  power, 
to  verify  whether  both  cases  generate  similar  requirements  or  not. 
The  methane  conversion  fraction  in  the  reformer  and  the  allowable 
maximum  solids  temperature  are  handled  as  discrete  decision 
variables,  for  further  structural  analysis  [27,28],  as  these  are 
believed  to  significantly  affect  the  thermomechanical  behaviour  of 
the  SRU.  The  lifetime  predictions  during  operation  at  constant 
system  power  are  compared  with  those  at  fixed  current  density, 
voltage  or  stack  power,  which  are  more  common  in  stack  testing. 

2.  Modelling  approach 

This  study  uses  a  one-dimensional  modular  model  of  an  SRU 
implemented  in  gPROMS  [29],  an  equation-oriented  process 
modelling  tool.  Thermal  boundary  conditions  are  imposed  to 
represent  the  conditions  in  an  SRU  embedded  in  a  stack.  The 
detailed  description  is  provided  in  Ref.  [25].  The  framework  is  the 
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Fig.  1.  (a)  View  of  the  FlameSOFC  SRU  design  and  one-dimensional  simplification  for 
the  calculations,  (b)  Description  of  domains  and  boundaries  in  the  model  [25].  7 Z: 
domains,  T:  interfacial  boundaries,  dU:  boundaries,  GDL:  gas  diffusion  layer,  MIC: 
metallic  interconnect,  area:  active  area,  inlet/outlet:  gas  manifold  zones. 
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design  developed  within  the  FP6  European  FlameSOFC  project.  Its 
objective  was  the  development  of  a  2.5  kWe  SOFC  stack  for 
combined  heat  and  power  (CHP)  application.  In  the  project,  fuel 
processing  was  based  on  thermal  partial  oxidation.  Here,  steam — 
methane  reforming  is  considered  in  addition,  because  of  the  higher 
achievable  efficiencies.  Fig.  1  depicts  the  actual  geometry,  along 
with  the  computational  domains  of  the  one-dimensional  model. 

The  electrochemical  model  considers  both  ohmic  and  non- 
ohmic  losses.  The  description  and  calibration  of  the  parameters 
on  experimental  data  from  two  different  segmented-cell  tests  was 
provided  elsewhere  [25].  The  main  features  are  recalled  hereafter. 

•  Ohmic  losses  comprise  the  ionic  resistivity  of  the  electrolyte, 
corrected  for  constriction  effects;  the  electronic  resistivity  of 
the  MIC,  its  oxide  layer  and  contact  resistance;  and  a  limited 
electronic  conductivity  of  the  electrolyte  that  induces  a  small 
leakage  current. 

•  On  the  cathode  side,  a  specific  modelling  approach  is  needed  for 
either  an  LSM-YSZ  or  else  lanthanum  strontium  cobaltite  ferrite 
(LSCF)  cathode.  In  this  study,  the  focus  is  on  LSM-YSZ  cathode.  A 
model  of  composite  electrode  solves  the  charge  balance  along 
with  mass  transport.  The  transfer  current  is  computed  following 
the  sequence  of  elementary  processes  proposed  by  van  Heuveln 
et  al.  [30],  adapted  for  a  low  coverage  of  adsorbed  oxygen  species. 
The  performance  of  the  LSM-YSZ  cathode  on  which  the  elec¬ 
trochemical  model  had  been  previously  calibrated  [25],  is 
significantly  lower  than  today’s  state-of-the-art.  For  instance,  the 
cell  potential  at  1073  K  and  0.3-0.7  A  cm-2  is  approximately 
10-50  mV  lower  than  that  reported  by  de  Haart  et  al.  [18].  This 
results  in  efficiencies  lower  than  expected  and  shorter  lifetime, 
but  does  not  preclude  the  analysis.  The  case  of  LSCF  cathode  has 
been  considered  for  comparison.  Here,  the  modelling  approach 
follows  that  proposed  by  Adler  et  al.  [31  ].  Experimental  data  from 
the  literature  are  used  to  compute  the  dependence  of  the  equi¬ 
librium  vacancy  concentration  on  temperature  and  oxygen 
partial  pressure  and  surface  exchange  properties  at  the  pore  wall. 
Direct  exchange  of  vacancies  with  the  electrolyte  is  assumed  [32]. 

•  On  the  anode  side,  the  dusty-gas  model  is  solved  in  one- 
dimension,  through  the  thickness  of  the  support,  along  with 
the  equation  of  continuity.  The  steam-methane  reforming  rate 
is  computed  by  the  kinetic  approach  of  Achenbach  and  Rien- 
sche  [33]  and  the  water-gas  shift  reaction  is  assumed  at 
equilibrium  in  the  anode.  Only  hydrogen  is  electrochemically 
converted  at  the  interface,  following  the  set  of  elementary 
processes  proposed  by  Zhu  et  al.  [34]. 

This  electrochemical  model  based  on  physical  principles  is 
compatible  with  the  implementation  of  degradation  phenomena. 
The  description  of  the  modelling  approach,  the  value  of  the 
parameters  and  the  analysis  of  the  local  behaviour  can  be  found  in 
Ref.  [20].  The  selected  degradation  phenomena  are  the  (i)  decrease 
in  ionic  conductivity  of  8YSZ,  (ii)  MIC  corrosion,  (iii)  anode  nickel 
particle  growth  and,  (iv)  chromium  contamination  and  (v)  forma¬ 
tion  of  zirconates  in  the  LSM-YSZ  cathode.  The  main  features  of  the 
degradation  models  are  summarised  as  follows: 

•  A  classical  percolation  model  predicts  the  reduction  of  the 
triple-phase  boundary  (TPB)  due  to  the  growth  of  nickel 
particles  in  the  anode.  The  semi-empirical  relation  for  the 
evolution  of  the  nickel  particle  radius  depends  on  temperature, 
steam  and  hydrogen  partial  pressure  and  tends  to  a  maximum 
value  due  to  the  mechanical  constraint  imposed  by  the  YSZ 
network.  In  the  model,  the  time  to  reach  the  plateau  value 
spans  from  1000  h  to  10,000  h  [20]  under  the  local  conditions 
found  in  an  intermediate  SOFC  stack. 


•  An  empirical  model  reproduces  the  experimental  data  on  the 
decrease  of  the  ionic  conductivity  of  8YSZ  samples  during  ageing 
in  air.  The  temperature-dependent  plateau  value  is  reached  in  the 
model  after  500  h— 2500  h  [20].  The  trend  in  the  evolution  after 
ageing  during  more  than  5000  h  is  not  included  in  the  data.  Errors 
may  arise  from  extrapolation.  This  degradation  process  affects 
both  the  electrolyte  and  the  8YSZ  phase  in  the  LSM-YSZ  cathode. 
Recovery  of  the  ionic  conductivity  is  not  allowed. 

•  The  maximum  SRU  temperature  in  intermediate-temperature 
SOFC  is  limited,  among  others,  by  the  oxidation  resistance 
and  high-temperature  mechanical  properties  of  the  metallic 
components.  This  is  not  included  in  the  present  study,  because 
a  consistent  set  of  data  on  the  temperature  dependence  of  the 
evolution  of  the  area  specific  resistance  (ASR)  of  coated  MIC 
could  not  be  found.  The  relation  by  Liu  et  al.  [35]  for  the  growth 
of  the  oxide  scale  under  a  Mn1.5C01.5O4  coating  that  is  used  in 
the  present  study  is  valid  at  1073  K  exclusively. 

•  The  modelling  of  chromium  contamination  of  LSM-YSZ  cathode 
assumes  the  progressive  blocking  of  the  active  sites  by  the  elec¬ 
trochemical  deposition  of  Cr203  from  Cr02(0H)2  (g)  released  by 
the  MIC  and  upstream  system  metallic  components.  This  process 
depends  upon  the  steam  and  Cr02(0H)2  (g)  partial  pressures  and 
assigned  to  the  local  overpotential  within  the  composite  cathode. 
Many  other  volatile  species  can  be  involved  [36-39].  For 
consistency,  a  small  amount  of  steam  is  needed  in  the  air.  The 
value  of  0.1%  used  here  is  typical  for  compressed  air,  commonly 
used  in  stack  tests  [40].  It  should  be  emphasised  that  depending 
on  the  humidity,  other  phenomena,  not  necessarily  related  to 
chromium  contamination,  may  predominantly  cause  the  degra¬ 
dation  of  the  LSM-YSZ  cathode  [41  ].  Chromium  contamination 
and  air  humidity  affect  LSCF  and  lanthanum  strontium  ferrite 
cathodes  as  well,  but  the  degradation  mechanisms  differ  [41—43]. 
This  restricts  the  use  of  the  chromium  contamination  modelling 
approach  developed  here  to  LSM-YSZ. 

•  Indicators  are  used  for  the  risk  of  reoxidation  of  the  anode  [21  ] 
and  formation  of  undesirable  zirconate  phases  (LZO,  SZO)  in 
the  LSM-YSZ  cathode.  For  the  latter,  the  thermodynamic  data 
on  the  critical  oxygen  partial  pressure  from  Liu  et  al.  [12]  is 
interpolated  and  compared  with  the  local  one,  computed  from 
the  local  overpotential  and  oxygen  gas  phase  concentration. 

Degradation  analysis  under  practical  operating  conditions 
requires  the  calculation  of  the  transport  of  volatile  chromium 
species  on  the  air  side  and  SOFC  system  considerations. 

2.1.  System  considerations 

The  aim  of  the  implementation  of  system  considerations  is  to 
capture  the  essential  implications  of  system  integration  on  the 
operation  of  the  SOFC  stack,  while  keeping  the  additional  compu¬ 
tational  load  as  low  as  possible  and  avoiding  the  complexity  of  an 


Fig.  2.  Schematic  view  of  a  typical  layout  of  SOFC  system  for  CHP  application. 
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Table  1 

Thickness  of  the  MEA  and  SRU  components. 


Electrolyte 

Compatibility  layer3 

Cathode 

Anode 

GDLa 

GDLC 

MIC 

Material 

8YSZ 

20GDC 

LSCF 

LSM-YSZ 

Ni-8YSZ 

Ni  foam 

Crofer22APUb  foam 

Crofer22APUb 

h  (m) 

7e-6 

7e-6 

60e-6 

540e-6 

le-3 

2e-3 

2e-3 

a  Implemented  in  the  case  of  LSCF  cathode. 
b  Commercial  denomination,  of  ferritic  alloy  [62]. 


extensive  system  integration  study.  Fig.  2  shows  a  schematic  view 
of  a  system  layout  for  CHP  application.  Previous  analyses  using  this 
system  configuration  [44]  have  demonstrated  a  large  system 
operating  envelope  with  sufficient  heat  integration  for  air  pre¬ 
heating  and  variation  of  the  methane  conversion  fraction  within 
the  reformer  in  the  range  of  conditions  considered  here.  Assuming 
the  system  is  well  thermally  integrated  as  has  been  shown  in  prior 
work  [44],  the  dominant  system  auxiliary  losses  that  dictate  system 
performance  are  those  due  to  the  air  blower  and  inverter. 

The  power  consumption  of  the  air  blower  per  SRU  in  the  stack  is 
computed  using  isentropic  relations  [44].  The  isentropic  efficiency 
is  0.8. 


Pblwr 


1 

blwr 


T 

T-  1 


RT, 


amb 


(1) 


The  calculation  of  the  stack  pressure  drop  on  the  air  side  is  based 
on  empirical  relations  for  the  permeability  of  metallic  foams  [45]: 


B0  =  81<r15exp(16.39n)  (2) 

A  simple  linear  approximation  is  used  for  the  dependence  on 
temperature  of  the  air  viscosity,  between  950  K  and  1150  K  [46]. 
The  variation  of  composition  along  the  channel  is  neglected, 
because  of  the  high  air  ratio: 


Qy  —  g  ^air  ifr  7^65  Over  7^iniet,  Parea?  ^outlet  (3) 

As  a  study  case,  the  pressure  drop  in  the  heat  exchanger 
network  in  Fig.  2,  needed  for  air  pre-heating  and  fuel  processing,  is 
assumed  to  be  twice  that  of  the  stack.  The  definition  of  the  specific 
system  power  is  based  on  the  SRU  active  area,  a  constant  inverter 
efficiency  of  93%  and  methane  fed  in  the  system. 
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Table  2 

Objectives,  values  of  the  discrete  decision  variables  and  range  of  the  continuous 
decision  variables  in  the  optimisation  problems. 


Objective 

System  efficiency  (eff) 

System  specific  power  (T^sys)3  (W  cm-2) 

Time  (h) 

Discrete  decision  variables 
Max.  SRU  temperature  (I<) 

Reformer  conversion  fraction 
Continuous  decision  variables5 
Fuel  utilisation 
Air  inlet  temperature  (K) 

a  Second  objective;  sequentially  varied  (see  text). 
b  Interval  for  optimisation. 


0.18;  0.35 
2000 

1100; 1125; 1150 
0.25;  0.50;  0.75;  0.99 


0.6-0.95 

600-1080 


2.2.  Transport  of  chromium  volatile  species 

At  the  SRU  level,  a  gas-tight  sealing  solution  yields  a  uniform 
steam  partial  pressure  on  the  cathode  side.  In  contrast,  the  content 
of  Cr02(0FI)2  released  by  the  MIC  is  governed  by  mass  transport 
limitations,  as  suggested  by  Opila  et  al.  [47].  Except  for  the  depo¬ 
sition  process  (see  Ref.  [20]),  the  present  approach  follows  that 
proposed  by  Wuillemin  [22],  the  main  features  of  which  are 
recalled  here. 

The  partial  pressure  of  Cr02(0H)2  over  the  oxide  scale  of  the  MIC 
is  assumed  at  equilibrium  and  computed  from  thermodynamic 
data  for  Cr203  [48,49]: 

log^ao2(OH)2  =  -2978-Tf1 -1.96  (6) 


*Cr02(0H)2  —  ^Cr02(0H)2  *XH20'x02  (^) 

In  the  fully  developed  region  of  a  laminar  flow,  the  analogy 
between  heat  and  mass  transport,  i.e.  Nu  Prnf  =  Sh  Scnt  [46], 
simplifies  the  one-dimensional  mass  conservation  equation  of 
Cr02(0H)2: 


„  Patm^Xc r02(0H)2  _  ^ vap  ,  st>deP 

airRTg  ay  “  JCCr02(OH)2  +  JtCr02(0H)2 


(8) 


JtCr02(0H)2 


Patm  NuDCv02(0H)2 
RT  L 


yeq 

ACr02(0H)2 


i  vl'z  vl'z 

dep  _  1  l°vcath *Cr02(0H)2*H20 

Cr02(OH)2  -  ha.r  An_F 

Z  =  heath 

J  Vcath (z)Atpb (z)2sin  h  Q d z 

z  =  0 


(9) 


(10) 


As  a  first  approximation,  Ucr02(0H)2  =  2.25  10-5  m2  s-1  [22]. 
This  approach  significantly  simplifies  the  real  situation  in  an 
SOFC  stack: 


•  In  ferritic  alloys  for  SOFC  application,  a  (Mn,Cr)304  spinel 
develops  on  top  of  the  oxide  scale  [50].  It  alleviates  the  release 
of  chromium  species,  due  to  the  thermodynamic  equilibrium, 
as  shown  by  computations  for  MnG^Oq  in  humid  air  [51  ].  The 
situation  is  further  complicated  in  a  coated  MIC.  The  use  of 
thermodynamic  data  for  the  equilibrium  partial  pressure  of 
Cr02(0FI)2  over  G2O3  therefore  yields  conservative  lifetime 
predictions,  even  though  the  scatter  in  the  thermodynamic 
data  induces  imprecision  in  the  range  of  one  order  of  magni¬ 
tude  [22,47,51]. 

•  Resistance  against  mass  transport  of  Cr02(0H)2  from  the  air 
channel  to  the  active  sites  where  it  is  electrochemically  depos¬ 
ited,  is  neglected.  This  assumption  is  conservative,  but  prevents 
the  use  of  the  model  for  reliable  guidance  in  the  choice  of  the 
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cathode  thickness  or  microstructure  or  current  collector  mate¬ 
rials  that  may  act  as  a  contamination-protecting  layer. 

The  present  modelling  approach  partially  captures  the 
complexity  of  SOFC  degradation,  because  the  underlying  physics  is 
highly  dependent  on  the  materials  and  conditions  and  the  under¬ 
standing  is  in  progress.  It  however  reproduces  the  correct  trends, 
without  any  specific  adjustment  [20],  but  cannot  comply  with  all 
the  observations  from  post-mortem  analyses  [18,52,53]  under  any 
tested  conditions.  For  such  an  achievement,  the  number  of  imple¬ 
mented  degradation  mechanisms  and  the  level  of  refinement  in 
their  description  must  be  increased,  which  requires  in  future 
a  more  intricate  interrelationship  between  experiments  and 
modelling.  We  believe,  for  example,  that  for  the  relevance  of  the 
chromium  contamination  model,  the  most  critical  limitations  arise 
from  the  lack  of  knowledge  on  the  deposition  process,  and  insuf¬ 
ficient  data  to  determine,  even  imperfectly  and  empirically,  the 
temperature  dependence  of  the  process. 


3.  Investigated  cases 

Apart  from  costs  and  application-related  considerations,  the 
end-user  specifies  the  electrical  power  output  of  the  SOFC  system 
that  must  be  delivered  at  the  highest  efficiency  and  power  density. 
The  target  electrical  power  of  a  given  SOFC  device  can  be  achieved 
by  manipulating  various  sets  of  operating  parameters.  A  sensitivity 
analysis  on  the  influence  of  each  single  parameter  on  the  initial  and 
long-term  performance  generates  a  large  amount  of  data,  the 
analysis  of  which  may  not  clearly  highlight  the  key  dependences. 
The  formulation  of  an  optimisation  problem  can  reduce  the 
complexity  and  amount  of  simulations  needed  for  a  useful  analysis. 

This  study  uses  the  simple  sequential  approach  for  multi-objective 
optimisation  applied  in  the  SOFC  field  by  Larrain  et  al.  [54].  The 
procedure  consists  in  performing  optimisations  on  one  objective,  the 
other  being  subjected  to  a  constraint  that  is  sequentially  modified  to 
cover  the  desired  range  of  variation.  A  second  run,  where  the  effective 
and  constrained  objectives  are  switched,  should  be  performed  for 
improved  reliability.  Computing  time  is  spared  here  for  extended 
variations  in  the  decision  variables,  rather  than  for  a  second  optimi¬ 
sation  run,  since  the  focus  is  on  the  identification  of  trends,  not  on  the 
effectiveness  of  the  multi-objective  optimisation  method. 

Optimisations  with  the  system  specific  power  and  efficiency  as 
conflicting  objectives  are  carried  out  separately  for  the  co-flow  (CO) 
and  counter-flow  (COU)  configurations. 

The  case  of  highest  system  efficiency  at  start  is  first  investigated. 
An  ideal  multi-objective  optimisation  includes  all  decision  vari¬ 
ables  to  provide  directly  design  guidelines.  In  reality,  both  thermo¬ 
electrochemical  and  structural  issues  govern  the  reliability  and 
durability  of  SOFC  devices.  Thermo-electrochemical  aspects  are 
included  in  the  one-dimensional  model,  to  the  limit  of  the  refine¬ 
ment  of  the  model  (see  Section  2  and  Refs.  [20,25]).  The  lack  of 
knowledge  on  mechanical  failures  in  SOFC  stacks,  and  the  variety  of 
their  nature,  prevent  their  simplification  to  indicators,  hence 
simple  implementation  in  the  optimisation  problem.  For  future 
structural  analysis  [27,28],  discrete  values  of  the  decision  variables 
that  mainly  govern  the  temperature  profile  in  an  SRU  embedded  in 
a  stack  are  enforced,  within  the  typical  range  for  intermediate- 
temperature,  anode-supported  SOFCs: 

•  The  maximum  allowable  SRU  solids  temperature  range 
considered  here  (1100  K,  1125  I<  or  1150  I<)  influences  the 
relative  contribution  of  residual  and  thermal  stresses  in  the 
different  components  of  an  SRU,  as  well  as  the  creep  and  rate- 
independent  plastic  behaviour  of  the  materials  [55]. 


•  The  variation  of  the  methane  conversion  fraction  in  the 
reformer  (25%,  50%,  75%  or  99%)  affects  the  magnitude,  position 
and  extent  of  zones  of  localised  higher  temperature  that  induce 
detrimental  stresses  in  the  SRU  components.  The  steam-to- 
carbon  ratio  is  two  for  all  cases. 

The  control  of  the  maximum  SRU  temperature  assumes  a  proper 
instrumentation  (direct  or  not)  of  the  stack  which  is  needed  if 
spatial  temperature  control  is  foreseen  [19].  The  air  inlet  temper¬ 
ature  determines  the  air  ratio  and  influences  the  temperature 
difference  over  the  SRU.  On  the  fuel  side,  the  fuel  utilisation  sets  the 
fuel  flow.  The  fuel  inlet  temperature  is  maintained  at  973  I<  in  all 
cases.  Table  1  lists  the  thicknesses  of  the  components  that  are 
identical  to  those  in  Ref.  [20].  Table  2  summarises  the  objectives 
and  decision  variables. 


System  specific  power  [W  cm  2] 

Fig.  3.  (a)  Optimised  initial  system  efficiency  (Eq.  (5))  and  (b)  SRU  potential  vs.  specific 
system  power,  for  the  co-flow  (dashed  lines)  or  counter-flow  (solid  lines)  configura¬ 
tions,  with  a  fixed  maximum  allowable  cell  temperature  of  1100  K,  1125  K  or  1150  K 
(black  to  light  grey),  and  methane  reforming  fraction  in  the  reformer  (PR)  of  0.25  (thick 
lines)  or  0.99  (thin  lines).  Optimisation  results  for  LSCF  cathode  are  provided  for 
comparison  (markers,  1100  K,  1150  K,  PR  =  0.25  and  counter-flow). 
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For  the  case  of  the  highest  system  efficiency  at  long-term, 
a  dynamic  optimisation  procedure,  where  the  time  dependence  ari¬ 
ses  from  the  degradation  phenomena  in  the  model,  is  performed 
using  the  same  objectives  and  decision  variables  as  for  the  steady- 
state  optimisation  procedure  (Table  2).  The  value  of  all  decision 
variables,  equivalently  control  variables,  is  kept  constant  during  the 
whole  ageing  period.  Indeed,  in  the  absence  of  a  refined  control, 
a  fixed  maximum  SRU  temperature  and  air  inlet  temperature  prevent 
a  drastic  modification  of  the  temperature  profile,  whereas  a  fixed  fuel 
utilisation  limits  the  risk  of  anode  reoxidation.  The  optimisations  are 
not  subjected  to  end  constraints  on  the  risk  of  formation  of  LZO/SZO 
and  anode  reoxidation.  These  indicators  are  monitored  for  analysis. 

The  time  horizon  influences  which  degradation  phenomena  are 
predominantly  alleviated  by  the  dynamic  optimisations.  A  lifetime 
criterion  is  not  straightforward  to  select  and  to  implement  in  the 
optimisation  problem.  A  fixed  time  horizon  represents  an  owner  life¬ 
time  requirement.  In  the  present  conditions,  a  value  of  2000  h  is  short 
enough  to  avoid  numerical  failures  in  the  high  system  specific  power 
regime,  yet  sufficient  to  capture  longer-term  characteristics,  since  the 
degradation  of  the  ionic  conductivity  of  8YSZ  and  of  the  anode  due  to 
nickel  particle  growth  are  close  to  their  plateau  at  this  time  [20].  Fixed 
gas  channel  geometry,  rather  than  constant  system  pressure  is 
enforced,  as  the  focus  is  on  the  operating  conditions  rather  than  on  SRU 
design,  and  for  consistency  with  structural  analysis  [27,28]. 

The  effect  of  operation  type  and  system  pressure  drop  on  the 
degradation  is  investigated  for  the  operating  conditions  for  the 


highest  electrical  efficiency  at  long-term  and  system  specific  power 
of  0.25  W  cm-2.  Three  scenarios  are  compared:  (i)  constant  current, 
(ii)  constant  SRU  voltage  and  (iii)  constant  system  specific  power. 
The  third  operation  mode  results  in  the  highest  current  levels  as 
degradation  lowers  the  voltage.  This  accelerates  degradation, 
predominantly  that  of  the  cathode  due  to  Cr  contamination  (see 
Section  4.2),  faster  than  (i)  and  (ii). 

Multi-objective  optimisation  results  commonly  exhibit  some 
scatter,  when  plotting  the  optimised  decision  variables  against  one  of 
the  objectives  [56,57].  For  clarity,  the  figures  that  provide  the  optimal 
values  of  the  decision  variables,  temperature  difference  and  SRU 
potential,  as  a  function  of  the  system  specific  power  use  trend  lines 
instead  of  actual  data.  The  maximum  error  is  provided  for  indication 
and  could  be  likely  reduced  by  a  second  optimisation  pass. 

4.  Results  and  discussion 

4.1.  Optimisation  of  the  system  electrical  efficiency  at  start 

Fig.  3  depicts  the  trade-off  between  system  efficiency  and  system 
specific  power  at  start,  and  corresponding  SRU  potentials.  The 
difference  between  the  trend  line  and  the  actual  data  is  in  the  worst 
case  less  than  1  %  and  3%  for  the  efficiency  and  potential,  respectively. 
Fig.  4  depicts  the  optimal  values  of  the  continuous  decision  variables 
and  the  temperature  difference  in  the  SRU.  The  efficiencies  fall 
within  the  range  reported  for  micro-CHP  devices  running  on  natural 


0.35 


System  specific  power  [W  cm  ] 


System  specific  power  [W  cm  2] 


Fig.  4.  Operating  conditions  to  achieve  the  highest  system  efficiency  at  start  depicted  in  Fig.  3:  (a)  air  ratio,  (b)  stack  air  inlet  temperature,  (c)  fuel  utilisation,  and  (d)  resulting 
temperature  difference  over  the  SRU,  for  the  co-flow  (dashed  lines)  and  counter-flow  configuration  (solid  lines),  with  a  fixed  maximum  allowable  cell  temperature  of  1100  K  (black) 
and  1150  K  (light  grey),  and  methane  conversion  fractions  in  the  reformer  of  0.25  (thick  lines)  and  0.99  (thin  line).  The  largest  discrepancy  between  the  trend  line  and  actual  results 
is  13%  (air  ratio). 
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gas.  The  specific  power  of  the  SRU,  i.e.  without  the  parasitic  power 
consumptions  of  the  air  blower  and  inverter,  and  the  average  current 
density,  ranges  from  0.22  to  0.41  W  cm-2  and  0.29  to  0.64  A  cm-2, 
respectively.  Because  of  the  conservative  performance  of  the 
considered  LSM-YSZ  cathode,  the  SRU  potential  is  low  for  this 
technological  solution.  For  comparison,  the  LSCF  cathode  yields 
a  significant  improvement  in  both  system  efficiency  and  SRU 
potential.  Among  the  discrete  decision  variables,  the  maximum  SRU 
solids  temperature  predominantly  influences  the  system  efficiency. 
Despite  the  low  SRU  potentials  and  high  fuel  utilisation,  anode 
reoxidation  affects  none  of  the  optimisation  results. 

The  gain  in  system  efficiency  due  to  a  decrease  of  the  system 
specific  power  is  less  significant  at  low  system  specific  power.  The 
display  is  restricted  to  a  value  of  0.2  W  cnrT2  because  of  the  effect  of 
the  leakage  current  in  the  electrolyte,  which  induces  a  decrease  of 
the  system  efficiency  starting  around  0.18  W  cm-2.  This  phenom¬ 
enon  has  been  implemented  because  in  the  data  used  for  the  model 
calibration,  open-circuit  voltages  slightly  lower  than  expected  were 
measured  [25].  This  has  been  previously  tentatively  ascribed  to  an 
increase  of  the  ionic  conductivity  of  the  thin  YSZ  electrolyte 
because  of  undesirable  doping  during  the  cell  production  [21  ].  Such 
phenomenon  has  not  been  observed  in  our  tests  of  more  recent 
cells,  for  which  sufficient  experimental  data  for  model  calibration  is 
however  currently  not  available.  The  simplified  modelling 
approach  of  this  phenomenon  hinders  a  detailed  and  reliable 
analysis.  Owing  to  the  high  electronic  ASR  of  the  electrolyte 
(50  Q  cm2),  the  leakage  current  density  is  lower  than  15  mA  cm-2  in 
most  cases,  which  results  in  a  decrease  of  the  system  efficiency  of 
1—2%  in  the  range  displayed  in  Fig.  3.  Compressive  gasket  materials 
are  not  completely  impervious  [58]  and  may  induce  a  similar  trend 
with  however  a  different  temperature  dependence. 

The  relationship  between  the  potential  and  specific  system 
power  follows  a  linear  trend  (Fig.  3).  In  contrast,  the  more 
pronounced  decrease  of  the  system  efficiency  at  high  specific  system 
power  is  due  to  that  of  the  fuel  utilisation,  which  must  be  lowered 
for  increased  current  density,  because  diffusion  limitations  issues  in 
the  thick  anode  support  become  more  significant.  The  relative 
parasitic  consumption  of  the  air  blower  further  amplifies  this  trend. 

4.1.1.  Discrete  decision  variables 

The  required  air  ratio  to  maintain  the  SRU  maximum  tempera¬ 
ture  within  specified  limits  is  high  (see  Fig.  4).  At  a  first  sight, 
internal  reforming  can  reduce  it  and  therefore  yield  a  benefit  on  the 
system  efficiency.  The  effect  of  methane  conversion  fraction  in  the 
reformer  depends  on  the  flow  configuration,  as  shown  in  Fig.  3.  It 
barely  influences  the  system  efficiency  in  co-flow,  whereas  internal 
reforming  (PR  =  0.25)  is  beneficial  in  counter-flow.  Fig.  5  qualita¬ 
tively  depicts  the  effect  of  methane  conversion  fraction  in  the 
reformer  on  the  temperature  and  current  density  profile  in  co-  and 
counter-flow  configuration.  In  co-flow,  the  location  of  the 
maximum  current  density  follows  a  trade-off  between  high 
temperature  and  lean  fuel  mixture.  At  high  fuel  utilisation,  it  is 
located  close  to  the  fuel  inlet  side  of  the  active  area,  which  is  the 
zone  cooled  by  the  endothermic  steam-methane  reforming  reac¬ 
tion.  Therefore,  the  air  inlet  temperature  is  comparatively 
increased,  as  shown  by  Fig.  4b,  to  avoid  a  very  low  electrochemical 
performance  in  this  area,  at  the  cost  of  a  higher  air  ratio,  hence  air 
blower  parasitic  consumption.  The  effect  of  internal  reforming  on 
the  system  efficiency  is  therefore  marginal  in  co-flow.  In  contrast,  it 
is  beneficial  for  the  system  efficiency  in  the  counter-flow  configu¬ 
ration,  because  the  cooling  effect  of  the  steam-methane  reforming 
reaction  acts  close  to  the  zone  of  highest  temperature,  i.e.  air  outlet 
side  of  the  active  area.  A  higher  air  inlet  temperature  is  therefore 
not  needed,  which  enables  a  significant  decrease  of  the  air  ratio, 
and  consequently  an  increase  of  the  system  efficiency. 


Fig.  5.  (a)  Temperature  profiles  over  the  SRU  and  (b)  current  density  profiles  over  the 
active  area  corresponding  to  operating  points  selected  from  Fig.  4:  maximum  cell 
temperature  of  1125  K,  system  specific  power  of  0.21  W  cm-2  and  0.29  W  cm-2  and 
methane  conversion  fraction  in  the  reformer  of  0.25  and  0.99.  j  in  (b)  refers  to  the 
average  current  density. 


4.1.2.  Continuous  decision  variables 

The  needed  air  ratio  scales  with  the  system  power  but  exhibits 
a  maximum  that  occurs  at  a  lower  specific  system  power  for  lower 
maximum  SRU  temperatures  (see  Fig.  4a).  At  low  specific  system 
power,  consequently  air  flow,  the  air  inlet  temperature  is  decreased 
because  the  in-plane  thermal  conductivity  of  the  SRU  is  kept 
constant.  Indeed,  Fig.  5a  shows  that  the  temperature  difference 
over  the  active  area  is  similar  at  0.21  W  cm-2  and  0.29  W  cm-2  in 
the  co-flow  configuration.  In  the  counter-flow  case,  only  part  of  the 
active  area  is  effectively  used  at  0.21  W  cnrr2.  The  temperature  at 
the  air  inlet  of  the  active  area  can  be  therefore  further  decreased,  as 
the  electrochemical  performance  in  this  zone  is  less  critical.  In  both 
flow  configurations,  the  lower  air  blower  consumption  ensuing  the 
reduction  of  the  air  ratio  increases  the  system  efficiency.  At  higher 
system  specific  power,  the  increase  of  the  air  inlet  temperature  is 
required  for  a  more  extensive  use  of  the  active  area,  but  is  balanced 
by  the  air  blower  parasitic  consumption. 

For  fixed  gas  inlet  temperatures  and  maximum  SRU  temperature, 
the  temperature  difference  over  the  SRU  scales  with  the  specific 
system  power,  because  of  air  flow  increases.  The  previous  discussion 
shows  that  this  is  not  the  case  for  practical  system  operating 
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conditions.  The  temperature  difference  is  the  highest  at  low  system 
specific  power,  because  of  the  low  corresponding  air  inlet  temper¬ 
ature.  In  co-flow  (counter-flow)  configuration,  a  complete  (low) 
methane  conversion  fraction  in  the  reformer  is  beneficial  in  terms  of 
system  electrical  efficiency,  but  yields  the  highest  temperature 
difference,  because  of  the  lower  air  inlet  temperature. 

The  dependence  of  the  optimised  fuel  utilisation  on  the  specific 
system  power  density  is  indirectly  affected  by  the  maximum 
temperature.  Diffusion  losses  exhibit  a  low  temperature  depen¬ 
dence  compared  with  other  contributions  to  the  losses  in  an  SOFC 
and  scale  with  the  current  density.  Higher  maximum  SRU 
temperature  yields  a  better  electrochemical  performance  of  the 
MEA,  hence  higher  voltage  and  lower  current  density  and  air  flow 
for  a  given  system  power.  The  fuel  utilisation,  consequently  system 
efficiency,  increases  as  shown  in  Fig.  4c. 


Fig.  6.  (a)  Typical  evolution  during  operation  at  constant  system  power,  of  the  system 
efficiency  (left,  Eq.  (5)),  SRU  potential  (left),  current  density  (left),  SRU  specific  power 
(left),  temperature  difference  over  the  SRU  (right)  and  relative  blower  power 
consumption,  as  a  fraction  of  stack  power  (right).  Counter-flow,  maximum  SRU 
temperature  of  1125  K,  methane  conversion  percentage  in  the  reformer  of  25%  and 
system  specific  power  of  0.25  W  cm-2,  (b)  Contributions  to  the  total  overpotential  in 
the  counter-flow  configuration,  at  a  constant  specific  system  power  of  0.29  W  cm-2, 
a  maximum  SRU  temperature  of  1125  K  and  methane  conversion  in  the  reformer  of 
25%.  The  Nernst  potential  is  computed  at  the  anode/electrolyte  interface. 


This  analysis  shows  that  comparing  flow  configurations  in  the 
light  of  identical  inlet  conditions  is  practically  of  limited  interest. 
Structural  aspects  may  limit  the  complete  use  of  the  benefits 
provided  by  operation  at  lower  specific  power,  higher  temperature 
and  adjustments  of  the  methane  conversion  fraction  in  the  reformer. 

4.2.  Long-term  system  electrical  efficiency 

The  less  stringent  operation  mode  for  an  SOFC  system  is 
constant  electrical  power  output  at  the  highest  efficiency.  Even 
though  this  best-case  scenario  excludes  any  cycling,  it  is  more 
demanding  than  most  stack  endurance  tests,  which  are  commonly 
performed  at  constant  current  density.  This,  along  with  the  rather 
low  performance  of  the  LSM-YSZ  cathode  used  for  the  calibration 
of  the  model  [25]  accounts  for  the  short  lifetimes  predicted  here. 

During  operation  at  constant  system  power  and  maximum 
allowable  temperature,  the  SRU  potential  decreases  because  of  the 


System  specific  power  [W  cm  2] 


System  specific  power  [W  cm  2] 


Fig.  7.  (a)  Initial  system  efficiency  (Eq.  (5))  and  (b)  cell  potential  vs.  system  specific 
power,  resulting  from  the  optimisation  of  the  system  efficiency  at  long-term  under 
operation  at  constant  system  power.  Co-flow  (dashed  lines)  or  counter-flow  (contin¬ 
uous  lines)  configurations,  with  discrete  maximum  allowable  cell  temperatures  of 
1100  K  (black)  and  1150  K  (grey),  and  methane  conversion  fraction  in  the  reformer  of 
0.25  (thick  lines)  and  0.99  (thin  lines). 
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degradation  processes  (see  Fig.  6a).  The  current  density  must  be 
ever  increased,  which  causes  higher  internal  polarisation  losses. 
The  air  ratio  increases  to  maintain  the  maximum  SRU  solids 
temperature  limit.  The  temperature  difference  increases  and  tends 
towards  the  value  determined  by  the  air  inlet  temperature  and 
maximum  allowable  temperature,  except  in  co-flow  configuration 
for  low  methane  conversion.  The  higher  air  blower  parasitic 
consumption  must  be  in  turn  compensated  by  an  increase  of  the 
current  density.  The  acceleration  of  the  regime  after  prolonged  use 
induces  the  end  of  life  of  the  device,  which  is  further  amplified  by 
the  sequential  activation  of  multiple  degradation  processes  (see 
Section  4.2.3).  As  discussed  in  detail  in  Ref.  [20],  the  degradation 
pattern  depicted  in  Fig.  6  has  been  observed  in  endurance  tests  of 
short  stacks  with  anode-supported  cells  based  on  an  LSM-YSZ 
cathode  [18].  Experimental  measurements  do  however  not  always 
exhibit  the  initial  decay  because  of  overlapping  activation 
processes,  which  are  not  implemented  in  the  model. 

Fig.  6b  depicts  the  contributions  to  the  total  overpotential  and 
their  evolution  upon  operation,  of  the  cathode,  anode,  electrolyte 
and  MIC.  The  contribution  of  the  cathode  dominates,  which  is  the 
case  for  all  the  investigated  cases  (not  depicted).  Among  the  other 
contributions,  that  of  the  MIC  exhibits  the  highest  relative  increase. 
Those  of  the  anode  and  electrolyte  are  comparable  (see  Part  II  [26]). 

4.2.1  Lifetime  improvement 

Fig.  7  depicts  the  initial  performance  that  results  from  optimi¬ 
sation  in  the  view  of  the  highest  system  efficiency  at  long-term.  In 


a  first  approximation,  the  optimiser  favours  the  SRU  potential  over 
the  initial  system  efficiency  (compare  Figs.  3  and  7).  The  effect  of 
the  maximum  allowable  temperature  and  methane  conversion 
fraction  in  the  reformer  is  similar  to  that  described  in  Section  4.1.  In 
co-flow  configuration,  a  complete  conversion  is  slightly  beneficial. 

Fig.  8  depicts  the  optimum  values  of  the  continuous  decision 
variables.  The  optimiser  reduces  the  fuel  utilisation  to  extend  the 
lifetime.  The  main  reason  is  not  the  degradation  of  the  anode,  as  in 
the  model  the  nickel  particle  radius  reaches  its  plateau  value  on 
most  of  the  active  area  before  2000  h  (see  Part  II  [26]),  but  the 
mitigation  of  the  progressive  increase  of  the  current  density  and  air 
ratio  during  operation,  while  allowing  a  slight  increase  of  the  air 
inlet  temperature  to  lower  the  cathode  overpotential.  This  results 
in  a  larger  difference  between  the  system  electrical  efficiency  at 
start  and  at  long  term  for  a  lower  maximum  SRU  temperature  and 
a  dependence  of  this  difference  on  the  second  objective,  the  system 
specific  power  (see  Fig.  7).  It  also  accounts  for  the  decrease  of  the 
temperature  difference  over  the  SRU  (compare  Figs.  4d  and  8d). 

Fig.  9  and  Table  3  provide  the  effect  of  the  operating  conditions 
and  their  optimisation  on  the  lifetime,  defined  here  as  the  time  to 
reach  a  drop  in  system  efficiency  of  20%.  The  computed  lifetime 
extension  is  provided  here  for  presentation  purpose.  The  model 
limitations  discussed  in  Section  2  and  incomplete  understanding  of 
the  degradation  phenomena  affecting  SOFCs  limit  their  applica¬ 
bility  to  all  cases.  The  counter-flow  configuration  along  with  the 
highest  maximum  SRU  temperature  (1150  K)  and  lowest  methane 
conversion  fraction  in  the  reformer  (PR  =  0.25)  results  in  the 


System  specific  power  [W  cm  2]  System  specific  power  [W  cm-2] 


Fig.  8.  Operating  conditions  to  achieve  the  highest  system  efficiency  at  long-term  as  a  function  of  the  system  specific  power  (Fig.  7):  (a)  air  ratio,  (b)  stack  air  inlet  temperature, 
(c)  fuel  utilisation,  and  (d)  resulting  temperature  difference  over  the  SRU,  with  fixed  maximum  allowable  cell  temperatures  of  1100  K  (black)  and  1150  K  (grey),  and  amount  of  fuel 
pre-reforming  (line  thickness).  The  largest  discrepancy  between  the  trend  line  and  the  actual  result  is  7%  (air  ratio). 
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System  specific  power  [W  cm"2] 

Fig.  9.  Effect  of  the  operating  conditions  depicted  in  Fig.  4  (markers)  and  Fig.  8  (lines) 
on  the  degradation  of  the  system  efficiency.  Operating  time  at  constant  system  power 
leading  to  a  decrease  in  efficiency  of  20%,  for  the  co-flow  (dashed  lines)  or  counter-flow 
(solid  line)  configuration,  with  a  maximum  allowable  cell  temperature  of  1100  K 
(black)  and  1150  K  (grey).  The  largest  discrepancy  between  the  trend  line  and  the 
actual  result  is  10%. 


highest  lifetime  and  averaged  system  efficiency  over  the  investi¬ 
gated  system  specific  power  range.  The  gain  in  lifetime  obtained  by 
manipulating  optimised  operating  conditions  are  listed  hereafter, 
from  the  most  to  the  less  significant  effect: 

•  System  specific  power:  factor  around  4  (0.21  W  cm-2  vs. 

0.28  W  cm'2). 

•  Maximum  SRU  temperature:  factor  around  2  (1150  I<  vs. 

1100  I<). 

•  Methane  conversion  fraction:  factor  around  1.3  (PR  =  0.25  vs. 

PR  =  0.99,  depending  on  flow-configuration). 

•  Flow  configuration:  factor  around  1.2  (counter  vs.  co-flow). 

In  the  conditions  treated  here,  the  degradation  of  the  cathode 
dominates.  Lower  cathode  overpotentials,  which  is  a  driving  force 
for  Cr  deposition  and  also  zirconate  formation  can  be  obtained  by 
adjusting  the  system  specific  power,  maximum  SRU  and  air  inlet 
temperature  and  flow  configuration.  The  benefit  of  decreasing  the 
overpotential  prevails  over  the  higher  vaporisation  rate  of  volatile 
chromium  species  from  the  MIC,  higher  relative  degradation  of  the 
conductivity  of  YSZ  and  higher  thermodynamic  risks  of  zirconate 
formation.  Manipulation  of  the  methane  conversion  fraction  is 
a  mean  to  lower  the  air  ratio,  consequently  the  increase  of  the 
system  parasitic  losses  to  maintain  the  specified  maximum  solid 


temperature  as  the  degradation  proceeds.  The  counter-flow 
configuration  limits  interactions  between  the  degradation 
phenomena  at  the  SRU  scale  compared  with  the  co-flow  configu¬ 
ration  (see  Section  4.2.3  for  anode  reoxidation  and  Part  II  for  anode 
and  electrolyte  degradation  promoted  by  the  evolution  of  the 
current  density  and  temperature  distribution).  The  aspects  are 
discussed  in  more  details  in  the  next  Sections. 

Our  structural  analyses  [27,28]  have  furthermore  shown  that  all 
these  manipulations  can  be  opposite  to  those  that  promote 
mechanical  reliability  and  durability,  resulting  in  trade-offs 
between  electrochemical  and  mechanical  considerations. 

4.2.2.  Effects  of  manipulations  of  the  operating  conditions  on  the 
lifetime 

Fig.  10  depicts  the  effect  of  methane  conversion  in  the  reformer, 
maximum  SRU  temperature,  assigned  system  power  and  flow 
configuration  on  the  evolution  of  the  system  efficiency.  The  pattern 
of  the  degradation  is  independent  of  the  operating  conditions.  The 
degradation  curves  of  the  system  efficiency  under  the  operating 
conditions  optimised  for  the  best  initial  performance  are  depicted 
for  a  specific  system  power  of  0.29  W  cm-2  (group  of  light  grey 
curves  A,  Fig.  10).  The  effects  of  the  maximum  SRU  temperature,  the 
methane  conversion  fraction  in  the  reformer  and  the  flow  config¬ 
uration,  are  consistent  with  Fig.  9  (see  Section  4.2.1). 

4.2.2 A.  System  specific  power.  Decreasing  the  specific  power  is 
expected  to  have  a  beneficial  effect  on  the  cost  of  electricity,  because 
it  promotes  higher  electrical  efficiency  and  the  gain  in  lifetime 
achieved  by  lowering  the  cathode  overpotential  largely  compen¬ 
sates  the  lower  power,  in  terms  of  produced  electrical  energy. 

Mechanical  reliability,  application-related,  manufacturing  and/ 
or  costs  considerations  that  are  not  included  in  the  present  opti¬ 
misation  problem  may  limit  the  advantages  of  these  adjustments. 
For  example,  operation  at  a  low  specific  power  to  extend  the  life¬ 
time  requires  to  stack  more  elements  to  achieve  the  targeted 
system  power.  This  increases  (i)  the  costs,  (ii)  the  risks  of  cell 
cracking,  due  to  the  larger  volume  of  brittle  ceramic  material  sub¬ 
jected  to  thermal  stress,  and  (iii)  gas  distribution  issues  among  the 
stacked  SRU.  The  lack  of  reproducibility  of  the  mechanical  proper¬ 
ties  of  the  component  materials  and  dimensional  tolerances 
amplifies  all  these  aspects.  In  this  view,  Table  4  and  the  comparison 
between  lines  and  markers  in  Fig.  9  suggests  the  significant  benefit 
of  targeting  higher  long-term  efficiency  on  the  lifetime,  by  a  factor 
up  to  3.5,  by  optimising  the  air  inlet  temperature  and  fuel  uti¬ 
lisation  for  given  system  specific  power,  maximum  SRU  tempera¬ 
ture,  methane  conversion  fraction  in  the  reformer  and  flow 
configuration. 

4.2.22.  Operating  temperature.  A  higher  operating  temperature 
promotes  extended  lifetime  by  reducing  the  cathode  overpotential, 
which  is  a  driving  force  for  both  the  electrochemical  deposition  of 
Cr203  at  the  TPB  and  formation  of  LZO/SZO  (20),  and  to  a  lesser 


Table  3 

Lifetime  extension  factors  obtained  by  manipulating  operating  conditions,  from  Fig.  8  (highest  long-term  system  efficiency).  The  lifetime  is  determined  as  the  time  to  reach 
a  decrease  in  system  efficiency  of  20%.  Inside  parentheses  is  the  averaged  system  efficiency  over  the  lifetime  and  produced  electrical  energy  (kW  h  cm-2). 


1100  K 

1125  K 

1150  K 

Co-flow,  0.29  W  cm-2,  PR  =  0.25 

Co-flow,  0.29  W  cm-2,  PR  =  0.99 

Co-flow,  0.21  W  cm-2,  PR  =  0.99 

1.00a  (0.441,0.71) 

1.06  (0.440,  0.75) 

4.04  (0.535,  2.08) 

1.53  (0.475,  1.09) 

1.73  (0.474,  1.23) 

5.92  (0.542,  3.05) 

2.14(0.493,  1.52) 
2.43  (0.493,  1.73) 
8.02  (0.539,  4.13) 

Counter-flow,  0.29  W  cm-2,  PR  =  0.25 
Counter-flow,  0.29  W  cm-2,  PR  =  0.99 
Counter-flow,  0.21  W  cm-2,  PR  =  0.25 

1.18(0.446,  0.84) 

0.92  (0.421,0.65) 

4.67  (0.545,  2.40) 

2.00  (0.480,  1.42) 

1.69  (0.466,1.20) 

7.14  (0.553,  3.68) 

2.63  (0.506,  1.87) 
2.18  (0.495,  1.55) 
9.71  (0.562,  5.00) 

a  Reference  case  for  the  calculation  of  the  lifetime  extension  factor. 
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Fig.  10.  Effect  of  (a)  the  amount  of  pre-reformed  methane  and  (b)  maximum  allowable 
cell  temperature  (co-flow:  PR  =  0.99,  counter- flow:  PR  =  0.25)  on  the  degradation  of 
the  system  efficiency,  for  the  co-flow  or  counter-flow  configuration,  with  the  operating 
conditions  depicted  in  Fig.  8,  except  the  light  grey  curves  (top,  A,  see  text)  (Fig.  4).  The 
occurrence,  from  thermodynamic  consideration,  of  anode  reoxidation  and  possible 
formation  of  LZO/SZO  is  indicated  by  empty  markers.  Constant  specific  system  power 
of  0.21  W  cm-2  and  0.29  W  cm-2. 

extent  by  reducing  the  parasitic  air  blower  consumption.  This  trend 
is  balanced  by  an  accelerated  increase  of  the  MIC  ASR,  due  to 
a  faster  growth  of  the  oxide  scale,  even  in  the  presence  of 
a  protective  coating  (see  Section  2),  higher  internal  vaporisation  of 


Table  4 

Lifetime  extension  factors  from  the  optimisation  of  the  operating  conditions,  i.e. 
ratio  between  the  conditions  depicted  in  Fig.  8  (highest  long-term  system  efficiency) 
and  Fig.  4  (highest  system  efficiency  at  start),  for  given  system  specific  powers  and 
maximum  SRU  temperatures.  The  lifetime  is  defined  as  the  time  to  reach  a  decrease 
in  system  efficiency  of  20%. 


1100  K 

1125  K 

1150  K 

Co-flow,  0.29  W  cm-2,  PR  = 

0.25 

2.88 

1.44 

1.67 

Co-flow,  0.29  W  cm"2,  PR  = 

0.99 

2.26 

1.73 

1.43 

Co-flow,  0.21  W  cm-2,  PR  = 

0.99 

1.05 

1.20 

1.08 

Counter-flow,  0.29  W  cm-2, 

PR  = 

0.25 

1.81 

1.63 

1.28 

Counter-flow,  0.29  W  cm-2, 

PR  = 

0.99 

3.46 

1.63 

1.15 

Counter-flow,  0.21  W  cm-2, 

PR  = 

0.25 

1.11 

1.06 

1.05 

volatile  chromium  species  from  the  MIC,  higher  relative  degrada¬ 
tion  of  8YSZ,  and  thermodynamic  risk  of  zirconate  formation.  In  the 
conditions  treated  here,  the  benefits  of  operation  at  higher 
temperature  prevail  over  the  drawbacks.  Investigations  with  the 
present  modelling  approach,  but  different  operating  conditions, 
suggest  that  the  maximum  lifetime  value  is  around  1180  K  [59], 
which  qualitatively  matches  the  operating  temperature  of  stacks 
with  electrolyte-supported  cells  and  coated  MICs,  e.g.  Ref.  [60].  The 
air  inlet  temperature  has  similarly  opposite  effects  on  the  degra¬ 
dation  of  an  LSM-YSZ  cathode.  Its  increase  yields  a  higher  inlet 
Pcro2(OH)2  and  higher  air  blower  parasitic  consumption  but  helps  to 
avoid  high  local  cathode  overpotentials. 

422.3.  Methane  conversion  fraction  in  the  reformer.  The  influence 
of  methane  conversion  fraction  in  the  reformer  is  opposite  in  the 
co-flow  and  counter-flow  configurations  (Section  4.2.1).  High 
conversion  values  are  detrimental  (beneficial)  in  the  counter-flow 
(co-flow)  configuration.  Despite  the  discrete  nature  of  this  decision 
variable  in  the  optimisation  problem,  the  dependence  is  likely 
monotonic  over  the  considered  range  of  variation,  as  suggested 
by  Fig.  10. 

4.2.3.  Progressive  activation  of  degradation  processes 

The  impact  of  the  reoxidation  of  the  anode  and  formation  of 
LZO/SZO  in  the  LSM-YSZ  cathode  on  the  electrochemical  perfor¬ 
mance  is  neither  modelled  nor  included  as  end-constraints  in  the 
optimisation  problem.  The  markers  in  Fig.  10  indicate  the  time  at 
which  these  phenomena  are  thermodynamically  activated  and 
would  accelerate  the  degradation.  The  reoxidation  of  the  anode 
does  not  occur  in  counter-flow  configuration  under  the  conditions 
studied  here,  owing  to  the  location  of  the  maximum  local  current 
density  close  to  the  fuel  inlet  zone  of  the  active  area,  where  the 
hydrogen  concentration  and  temperature  are  the  highest.  In  co¬ 
flow,  it  is  activated  after  the  formation  of  zirconates  and  can  be 
practically  avoided  by  adjusting  the  fuel  utilisation.  It  arises  from 
the  progressive  increase  of  the  current  density  and  occurs  first  at 
the  anode/electrolyte  interface  in  the  zones  where  the  fuel  is 
depleted.  The  lifetime  defined  as  the  time  at  which  anode  reox¬ 
idation  may  take  place,  qualitatively  follows  the  same  trend  as  in 
Fig.  9  (not  depicted). 

Fig.  11  depicts  the  dependence  of  the  operating  time  at  which 
zirconate  can  form  in  the  operating  conditions.  The  effect  of  the 
operating  temperature  is  governed  by  a  trade-off  between  (i)  lower 
cathode  overpotential  and  (ii)  the  thermodynamic  promotion  of 
the  formation  of  zirconate  in  the  LSM-YSZ  system  at  high 
temperature,  and  oxygen  partial  pressure  at  the  cathode/electrolyte 
interface  (see  Part  II  [26]).  The  effect  of  (i)  dominates  that  of  (ii)  over 
the  investigated  range.  For  a  given  maximum  SRU  temperature  and 
system  specific  power,  the  co-flow  configuration  provides  better 
overall  resistance,  since  the  zone  of  high  cathode  overpotential  is 
subjected  to  lower  temperature.  Even  though  the  formation  of  LZO / 
SZO  is  not  included  in  the  optimisation  problem,  it  is  delayed  by  the 
reduction  of  the  overpotential  to  alleviate  the  chromium  contam¬ 
ination  of  the  cathode. 

The  aforementioned  results  highlight  the  importance  of  the 
cathode  overpotential  on  the  degradation  of  the  performance. 
Fig.  12  shows  the  maximum  local  cathode  overpotentials  induced 
by  the  operating  conditions  optimised  for  the  highest  system  effi¬ 
ciency  at  start  (Figs.  3  and  4)  and  long-term  (Figs.  7  and  8).  The 
trade-off  between  temperature  and  anode  gas  composition  in  the 
co-flow  case  induces  significantly  higher  cathode  overpotential 
than  in  the  counter-flow  case.  The  optimisation  of  the  long-term 
system  efficiency  clearly  favours  lower  cathode  overpotential,  by 
a  slight  increase  of  the  air  inlet  temperature  and  reduction  of  the 
fuel  utilisation  that  alleviates  the  constraint  of  the  distribution  of 
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System  specific  power  [W  cm"2] 

Fig.  11.  Operating  time  at  constant  system  power,  before  the  possible  formation  of 
zirconates  in  the  LSM-YSZ  cathode,  in  an  SRU  operated  in  co-flow  (dashed  lines)  or 
counter- flow  (solid  lines)  configuration.  Conditions  from  Fig.  8.  The  maximum  allow¬ 
able  cell  temperature  of  1100  K  and  1150  K  is  indicated  by  grey  scale  variation  and  the 
methane  conversion  fraction  in  the  reformer  by  line  thickness.  The  largest  discrepancy 
between  the  trend  line  and  the  actual  result  is  13%. 


System  specific  power  [W  cm  2] 


the  gas  composition  over  the  active  area  on  that  of  the  current 
density.  The  effects  of  the  maximum  SRU  temperature  and  methane 
conversion  fraction  in  the  reformer  are  minor  at  the  initial  time,  but 
are  amplified  by  ageing.  The  dependence  of  the  maximum  cathode 
overpotential  on  the  decision  variables  and  system  specific  power 
follow  similar  trends  than  the  lifetime. 

4.3.  Effect  of  system  operation  and  cathode  performance 

Endurance  tests  of  SOFC  stacks  are  commonly  performed  at 
constant  current  density  [18,61],  which  is  favourable  for  the  stack 
lifetime  but  does  not  fully  reflect  field  conditions.  Fig.  13  shows 
the  difference  between  operation  at  constant  potential,  current 
density  and  system  specific  power.  Testing  conditions  strongly 
influence  the  observed  degradation  and  evolution  of  the  air  ratio. 
The  reoxidation  of  the  anode  and  formation  of  LZO/SZO  do  not 
occur  during  operation  at  constant  potential,  owing  to  the  small 
variation  of  the  cathode  overpotential  and  monotonic  decrease  of 
the  current  density.  For  the  same  reasons,  the  air  ratio,  hence 
temperature  difference,  decreases.  The  reoxidation  of  the  anode 
may  occur  at  constant  current  density,  in  the  co-flow  configura¬ 
tion,  due  to  the  redistribution  of  the  current  density.  As  the 
cathode  performance  decreases,  the  zone  of  highest  current 
density  is  shifted  towards  zones  of  higher  temperature  and  lower 
hydrogen  molar  fraction.  This  is  not  the  case  in  the  counter-flow 
configuration. 


System  specific  power  [W  cm-2] 


Fig.  12.  Maximum  local  cathode  overpotential  at  the  initial  time  (a,  b)  and  after  3000  h  of  operation  at  a  constant  system  power  (c,  d),  under  the  conditions  depicted  in  Figs.  4  and  8, 
for  the  co-flow  (dashed  lines)  or  counter-flow  (solid  lines)  configuration,  with  different  discrete  maximum  allowable  cell  temperatures  (grayscale  variations)  and  methane 
conversion  fraction  in  the  reformer  (line  thickness).  The  largest  discrepancy  between  the  trend  line  and  the  actual  result  is  6.5%. 
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Fig.  13.  Comparison  of  (a)  the  evolution  of  the  SRU  potential,  current  density,  specific  power,  (b)  temperature  difference  and  air  ratio,  during  operation  at  constant  potential  (light 
grey),  current  density  (grey)  and  system  power  (black),  for  the  co-flow  (PR  =  0.99)  and  counter-flow  (PR  =  0.25)  configurations.  The  maximum  cell  temperature  is  1125  K.  An 
identical  initial  cell  potential  is  enforced  among  the  operating  scenario,  yielding  an  initial  system  specific  power  of  0.25  W  cm-2.  Empty  markers  indicate  the  occurrence  of  anode 
reoxidation  and  possible  formation  of  LZO/SZO. 


The  combined  influence  of  the  air  ratio,  air  inlet  temperature 
and  fuel  utilisation  identified  in  Section  4.1.1  underscores  the 
importance  of  including  system  considerations  in  lifetime  analyses. 
The  parasitic  consumption  of  the  air  blower  is  related  to  the  stack 
pressure  drop.  Fig.  14  shows  that  the  system  pressure  drop  has 
a  significant  effect  on  the  degradation  of  the  system  efficiency.  For 
this  comparison,  the  stack  specific  power,  rather  than  the  system 
specific  power  is  kept  constant  to  enforce  identical  initial  SRU 
conditions,  i.e.  overpotentials.  A  stack  and  interconnection  design 
for  extended  lifetime  must  therefore  consider  pressure  drop  in  the 
light  of  (i)  the  supply  of  reactant  to  the  active  sites  and  current 
constriction  effects  in  the  interconnection  system,  (ii)  the  need  of 
a  homogeneous  in-plane  gas  distribution,  (iii)  the  benefits  of  lower 
system  parasitic  consumptions  and  (iv)  the  mechanical  stability  of 
the  GDL  and  sealants. 


Fig.  14.  Effect  of  the  system  pressure  drop  on  the  evolution  of  the  system  efficiency, 
SRU  potential  and  blower  consumption  per  SRU,  during  operation  at  a  constant  system 
specific  power  of  approximately  0.25  W  cm-2,  the  identical  initial  SRU  power  is 
enforced  among  the  nominal  and  high  pressure  drop  cases.  Co-flow:  PR  =  0.99, 
counter-flow:  PR  =  0.25. 


5.  Conclusion 

The  developed  modelling  approach  aims  at  supporting  stack 
and  system  design  and  operation  for  extended  durability.  This 
study  has  placed  the  emphasis  on  the  effect  of  the  operation 
conditions  on  the  degradation  of  SOFC  stacks  integrated  in 
a  system,  under  practical  constant  system  electrical  power  output 
mode.  Optimisation  has  shed  light  on  critical  decision  variables, 
their  inter-dependences  and  the  differences  between  the  require¬ 
ments  for  the  highest  system  electrical  efficiency  at  start-up  and  at 
long-term. 

The  degradation  processes  detrimentally  interact  in  an  SRU. 
Under  practical  operating  conditions,  their  sequential  activation 
shortens  the  lifetime  of  the  device.  The  loss  of  performance  due  to 
chromium  contamination  of  an  LSM-YSZ  cathode  can  activate  the 
formation  of  LZO/SZO  insulating  phases.  Reoxidation  of  the  anode 
can  occur  in  the  co-flow  configuration  if  the  fuel  utilisation  is  not 
adjusted  upon  operation,  because  of  the  increase  of  the  current 
density  to  deliver  a  constant  system  electrical  power  output.  It  did 
not  occur  during  the  simulations  in  counter-flow  configuration.  The 
system,  through  the  air  blower  parasitic  consumption,  also  plays  an 
important  role  in  the  system  efficiency  degradation. 

Besides  material  improvements,  the  manipulation  of  the  oper- 
g1  ating  parameters  is  an  effective  and  practical  means  to  improve  the 
3  durability,  which,  like  system  electrical  efficiency,  scales  inversely 

§5  with  the  system  specific  power.  The  gain  in  lifetime  largely 

£  compensates  for  the  decrease  in  specific  power,  in  terms  of 
produced  electricity.  The  lower  maximum  cathode  overpotential  in 
counter-flow  and  the  effect  of  internal  reforming  on  lowering  the 
8-  parasitic  air  blower  consumption  is  beneficial  since  it  alleviates  in 
the  long-term  cathode  chromium  contamination,  risk  of  anode 
reoxidation  and  zirconate  formation  in  the  cathode.  In  reality, 
mechanical,  application,  manufacturing  and/or  costs  consider¬ 
ations  may  limit  the  advantages  of  these  adjustments.  Examples 
that  are  not  considered  in  the  model  are  gas  flow  distribution 
problems  arising  from  the  need  to  increase  the  number  of  SRU  in 
the  stack  and  larger  irreversible  deformation  generated  at  high 
temperature.  Future  stress  analysis  is  needed  to  better  characterise 
the  trade-off  between  electrochemical  and  mechanical  degrada¬ 
tion.  In  this  view,  the  separate  optimisation  of  the  air  inlet 
temperature  and  fuel  utilisation  in  the  view  of  long-term,  instead  of 
initial  system  electrical  efficiency,  further  extends  the  durability 
(factor  up  to  3.5). 
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The  present  study  only  partially  addresses  the  complexity  of 
SOFC  degradation,  because  the  underlying  physics  is  highly 
dependent  on  the  materials  and  application,  and  deeper  under¬ 
standing  is  in  progress.  It  however  reveals  the  potential  for 
considerable  durability  improvements.  To  extend  the  knowledge 
and  push  further  the  current  limitations  of  the  SOFC  technology, 

•  Relevant  analyses  must  include  the  interactions  between  the 
various  phenomena  that  affect  SOFC  degradation. 

•  Control  and  operation  strategies,  as  well  as  system  integration, 
must  be  integrated  with  lifetime  analysis. 

•  Modelling  and  experimental  tasks  must  be  carried  out  in 
synergy. 
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Nomenclature 

Latin  letters 

N  molar  flow  rate  (mol  s-1) 

V  specific  power  (W  m-2) 

V  volume  (m3) 

rate  of  production  of  species  i  (mol  m-3  s-1) 

LHV  lower  heating  value  (J  mol-1) 

A  area  (m2) 

ATpb  specific  area  at  the  triple  phase  boundary  (m2  m-3) 

B0  permeability  in  porous  medium  (m2) 

Df-  effective  bulk  diffusivity  of  binary  pair  in  porous  medium 
J  (m2  s-1) 

Ea  activation  energy  (J  mol-1) 

F  Faraday's  constant  96,485  (C  mol-1) 

h  thickness  (m) 

i0  exchange  current  density  (A  m-2) 

k0  kinetic  constant 

I<e q  equilibrium  constant,  dimensionless 

L  length  (m) 

mcrp  creep  stress  exponent 

n  porosity 

nt  exponent 

Nu  Nusselt  number 

P  power  (W) 

p  pressure  (Pa) 

patm  atmospheric  pressure  101,325  (Pa) 

R  universal  gas  constant  8.314  (J  mol-1  I<-1) 

Iamb  ambient  temperature,  298  (I<) 

u  gas  velocity  (m  s-1) 

xj  mole  fraction  of  species  i 

Greek  letters 
e  strain 

eeff  efficiency 

7 1  overpotential  (V) 

y  ratio  of  specific  heats  ( cp/cv ) 

pg  viscosity  (Pa  s) 

p  density  (kg  m-3) 

cr  stress  (Pa) 


Indices 

V  degradation 

blwr  blower 

cath  cathode 

eq  equivalent 

o  initial,  dense  state 

out  outlet 

s  index  for  solid  parts,  individual  components  or  averaged 

structure 
sys  system 

Superscripts 
*  foam 

crp  creep 

dep  deposition 

vap  vaporation 

Acronyms 

ASR  area  specific  resistance 

CFIP  combined  heat  and  power  generation 

COU  counter-flow 

CO  co-flow 

EIS  electrochemical  impedance  spectra 

GDL  gas  diffusion  layer 

LSCF  lanthanum  strontium  cobaltite  ferrite 

LSM  lanthanum  strontium  manganite 

LZO  lanthanum  zirconate 

MEA  membrane  electrode  assembly 

MIC  metallic  interconnect 

SRU  single  repeating  unit 

SZO  strontium  zirconate 

TPB  triple  phase  boundary 

YSZ  yttria-stabilised  zirconia 
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